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Abstract. Using Cholesky decomposition of the moment matrix of measure 
a, we define the set of polynomials orthonormal with respect to a. We give 
some properties of them including those related to Christoffcl function. We are 
also able to give simple formula for expansion of monomial x n in orthonormal 
polynomials and define general algorithm for obtaining the so called lineariza- 
tion coefficients. 

Considering two measures a and 5 and two sets of polynomials orthogonal 
with respect to them we are able to give general formula for the connection 
coefficients between the two sets of polynomials. Moreover if a << <5 and 
Radon-Nikodym derivative da / d& is square integrable with respect to dS then 
we expand da/d5 in Fourier series of polynomials orthonormal with respect to 
<5. We illustrate developed theory by providing yet another proof of the famous 
Poisson— Mehler expansion formula. 

We also find general expressions for the power coefficients of orthogonal 
polynomials and moments in terms of the coefficients of 3-term recurrence at 
least for symmetric distributions. 

In all these considerations the new interpretations of the inverse of the 
moment matrix are presented. 



1. Introduction 

Let us first make some remarks concerning notation, a, /3, . . . will denote positive 
measures. In order to be able to use sometimes probabilistic notation we will assume 
that all considered measures are normalized. Integrals of integrable function / with 
respect to measure say a will be denoted by either of the following notations 

f(x)da(s), J fda,Ef, Ef(X), E a f(X), 

depending on the context and the need to specify details. Above X denotes random 
variable with distribution a. Probability theory assures that it always exist. 

Matrices and vectors (always columns) will be generally denoted by bold type let- 
ters. The most important vector and matrix are X„ = (1, x, . . . , x n ) T (T— transposition) 
and 

(1.1) M n (a) = [^(q)]. h jn , 

where m n {a) = J x n da(x). Another words M„(a) = fiaX.n.Xj'. Let us also intro- 
duce also vectors consisting of successive moments 

= (l,...,m„). 
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Let us remark immediately that the matrix M« is the main sub-matrix of the matrix 
M n +i. Let us also introduce sequence 

(1.2) A„(a) = detM n (a), 

n > 1, of determinants of matrices M„(a). 

In order to avoid repetition of assumption we will assume that matrices M„ exist 
for all n > 0. Obviously (0,0) entry of the matrix M„ is equal to 1. 

We know that given measure a such that all moments exist one can define set 
of polynomials {p n {x, a)} n> _ 1 with p-i(x,a) = 0, po(x,a) = 1 and such that 

p n (x,a)pm{x,a)da(x) = S n<m , 

where 8 n:in denoted Kronecker's delta. Moreover if we declare that all leading coeffi- 
cients of polynomials p n (x, a) are positive then coefficients n n ,i(a) of the expansion 

n 

(1.3) p n (x, a) = 7r TM (a) x\ 

i=0 

are defined uniquely by the measure a. According to our convention, later we will 
drop dependence on a, if measure a is clearly specified. 

Let us define vectors P„(x) = (po(x), ■ ■ ■ ,p n (x)) T and the lower triangular ma- 
trix II n with entries Tti,j- Of course we set nij = for j > i. We obviously have: 

(1.4) p„(.T) = n„x„. 

To continue introduction of notation let X n ,i{a) denote coefficients in the following 
expansions: 

n 

(1.5) x n =^2\ n ,i(a)pi(x,a) . 

Consequently let us introduce lower triangular matrices A„ with entries A^.j if i > j 
and otherwise. 
We obviously have: 

(1.6) X„ = A„P„(x), n„A„ = A„n„ = I„, 

where I„ denotes (n + 1) x (n + 1) identity matrix. 

Since polynomials {p n } are orthonormal then there exist two number sequences 
{a„} , {b n } such that polynomials {p n } satisfy the following 3-term recurrence: 

(1.7) xp n {x) = a„ +1 p n+1 (x) + b 
with ao = and n > 0. We know also that 

(1.8) a n = 7Fn ~ 1 -" 1 ; h n — xp n (x)da(x) 7 

consequently that bo = mi. For details see e.g. p] or [3J. 

Combining (|1.7j) and (|1.8|) we get the following set of recursive equations to be 
satisfied by coefficients Tr n j ■ 

(1-9) Cln+l^n+lM + b n TT n fi + a n 7Tn~lfi = 0, 

(1-10) a n+1 TV n+1 j + b n TV n j + a n Tr n -l.j = 7T n ^_l, 

for n > 0, j — 1, . . . , n , remembering that 7r nj = for j > n. 
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Remark 1. As it follows from formula (2.1.6) of [3] coefficients Tr n ,i can be ex- 
pressed as determinants of certain submatrices built of moment matrix M„. In 
particular denoting by Dn'^ the determinant of a submatrix obtained by removing 
row number i + 1 and column number j + 1 of the matrix M„. We have Tr n ^ = 

(-l)"-^ n) /V A « A «-i- 

Let us also consider family of associated polynomials {q n (x)} n> _ 1 . As it follows 
say pQ or [7] they satisfy the same 3-term recurrence but with different initial values. 
Namely we assume that q~i(x) = — 1 and qo(x) — 0. Following (|1.7[) we see that 
then qi(x) — l/a\. One knows also (see e.g. [7]) that 

q n (x) = / da[y). 

J x-y 

Now since (x k - y k )/(x - y) — x k ~ x + x k ~ 2 y + . . . + y fe ~\ p„(x) = 2™ =0 n n ,jX j 
and J y k da (y) — we deduce that 

n j— 1 n— 1 n 

(1.11) q n (x) = ^7r„j ■^mj_i_fcX fe = ^x k ^2 'Knjmj-x-k- 

j = l k=0 k=0 j=k+l 

Let us define also (n + 1)— vector Q n (x) — (0, qi(x), . . . , q n (x)) T . 

It should be stressed that our results in fact mostly concern the so called " trun- 
cated moment problem" that is we in fact assume that we know finite number (say 
In + 1 including moment or order 0) of moments of some distribution. That is in 
many cases the assumption that the matrix M„ exists for all n will not be needed. 
Then we derive n polynomial {f>i}™ =0 that are mutually orthogonal, we find coef- 
ficients of expansion of x l in terms of these polynomials as well as we derive all 
so called linearization coefficients i.e. coefficients of the expansions pi(x)pj(x) in 
polynomials {pi}™ =0 for all i+j< n. 

Given two distributions (say a and ft) and 2 respective moment sequences we 
are able to derive all so called "connection coefficients" i.e. coefficients of the 
expansion of say pj (x,5) in {pi(x, a)} 3 i=0 and conversely. Due to very efficient 
numerical algorithms of Cholesky decomposition and inversion of lower triangular 
matrices all these calculations can be done within seconds using today's computers. 

Of course we present also results that require existence of all moments. These are 
some limit properties of arithmetic averages of orthogonal polynomials and more 
importantly results concerning expansions of Radon-Nikodym derivatives of one 
distribution with respect to the other (see (|2.8|l ). 

The paper is organized as follows. In the next Section [5] are our main results. 
As we think particularly interesting are the ones concerning connection coefficients 
presented in Subsection 12.11 containing not only formula for the connection coeffi- 
cients between two sets of orthogonal polynomials related to two measures but also 
expansion of Radon-Nikodym derivative of one measure with respect to the other 
in a Fourier series of orthogonal polynomials related to one of the measures. In- 
teresting seems also Subsection 12.21 presenting general formula for the linearization 
coefficients. Longer and uninteresting proofs are shifted to Section [3] 
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2. Cholesky decomposition and main results 

Our basic tool in what follows is the so called Cholesky decomposition of the 
symmetric, positive definite matrix. Below we collect some of the properties of 
Cholesky decomposition in the following simple proposition. 

Proposition 1 (Cholesky). Suppose positive normalize measure a has support of 
infinite cardinality and J x 2N da for some N > 0. Then 

i) V N > n > 1 there exists unique real, non-singular lower triangular matrix 
L n (a) such that M„(a) = L n (a)L„ (a). 

ii) entries of matrix L^r can be calculated recursively 



(2.1) 



\ 



n-1 fc-1 

m2n — ^ ln+l,k = { m n+k+l — ^ l n +l ,jh,j)/h,ki 

3=0 j=0 



with Iq.q = 1 for n = 0, . . . , N. Entries l n , n have also the following interpretation: 

(2.2) r A " 



n,n 



where sequence {A n } is defined by M.ty) . In particular we have: 



(2.:-!) = ^/m 2 -mf, h,2 = ^rrn- 

(tth+i - m^x) 



(m 3 - mim 2 ) 
(m 2 - mf) 



2 



(2.4) Z jj0 = mj, = 



/i. 



, n cN , 1 , (m 3 - m 2 m 1 m l+1 - m.mi 

(2.5) Z ij2 = (m i+2 - 171^2 t or ), 

«2,2 (m 2 -mf) 

i = l,...,N, 

Hi) V < i, j < N 

min(tj') 

m i+j = h,klj,k- 
k=0 

Proof, i) Follows existence and uniqueness of Cholesky decomposition (see e.g. The- 
orem 8.2.1 of [5]) and the fact that if support of positive measure is infinite then 
matrix M^r exists, is symmetric and positive definite. Besides by Cauchy Theorem 
we have A„ = (dcti„) 2 = n"=o'i,r Since A„_i = YljZo we S et our assertion, 
ii) Follows one of the algorithms of obtaining Cholesky decomposition (so called 
Cholesky-Banachiewicz algorithm). □ 

We have obvious observations: 

Proposition 2. Let M n and M" 1 , n > be respectively sequence of moment 
matrices and the sequence of its inverses of some measure a. Let us denote M^ 1 
= [/4" ]o<i,j<n i- e - that /if 1 ) is entry of the matrix M^ 1 . Let L„ be defined 

by it sequence of lower triangular matrices forming Cholesky decomposition, then 

i) Mn > 0, II„ = L" 1 , A„ = L„. That is A„A^ = M„ and II^II„ = M" 1 in 

particular ELmax(ij) ^k,i^k,j = 4d and EfcS) Kk\j,k = m i+j . 

ii) pT[x)P n {y) = E?=oPi(x)Pi(y) = X^M^ X Y W , thus X^M"^ is the re- 
producing kernel and l/X^M r ^ 1 X„ is the Christoffel function of the measure a. 
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Consequently 



1 + . . . + a; 2 " - X^M,T 1 X„ ~ 1 + . . . + x 2n ' 

where Ao. n arcd A„. n are respectively the smallest and the largest eigenvalues of the 
matrix M„ . 

Hi) J Pl(x)M n P n (x)da(x) = E"=o m 2i = A ,„ + . . . , A„,„, 

W ^/ 2T P„( e lt )P^( e - lt )* = M- 1 consequently, ^ E7=o it H^f ^ = 



MM- 1 ) = Ei=o 1 Ao 



) . 7 i. 1 



^ ^fi? ^ ^"=1 I* (°)! 2 = ^i<iJ<n ''<>< >'»./ : < ^S^ 1 , «^E"=i <&(0) Pj -(0) 

En (n) 



wzj , , s , 1 r; — rrrr V!? n Pi aO — ^ 0, a— a.s. as n 

y Vn+l log J (n+2) i-d=0-^'V ' ' ' 



00. 



Proof. Is shifted to Section |3l □ 

Remark 2. Assertion iv) was shown in [2]. We present it here for completeness. 

Remark 3. Notice that we have Ej=o Pj ( z ) — Ej=o t jL Oj z '' ■ Now since Ej>o biWT 
= /lIqq by assetion v) we see that in the case of determined moment problem it might 
happen that Ej>o |Pj(0)| 2 = 00 an d consequently that one of the Nevanllina func- 
tions (D(z) as defined by Theorem 4-9 in ^\) might not be determined in this case. 

Remark 4. Notice that from vi) it follows that if the problem is indeterminate i.e. 
when $^=i \lj (0)1 ^ 00 ( see either Theorem 2.17 or [6],) we get estimate of speed 
of the divergence of X n . n to infinity. 

As a corollary we have the following observations: 

Corollary 1. Coefficients a n and b n : n > defining 3-term recurrence are related 
to the moment matrix by the formulae: 

(0 R\ 2 _ A n A ra _2 _ An-i A„_ 2 

^Z.OJ d n — . 2 7 — a ^n+l,n^n,n » 'n,n-Pn-l,n-li 

/or n > 2 im'f/i ao = 0, a\ — A2 = m2 — m\. 

Proof. Following (|1.8I) and Proposition [5] i) we deduce a\ = tt 2 1 _ 1 n _ 1 /ir 2 l n . Since 
TTn.n = C,n we a PPly (|2.2I) . To get formula for 6„ first we observe that (i,i — 1) 

entry of the the inverse of the lower triangular matrix L„ = [lij]i=o n,i=o,...,< is 

equal to — ; ~ 1 , - ■ Besides dividing both sides of (| 1 . 10[) with j = n by 7r„ n we 
get: 

■ h £>„ = ' . 

^"+1." _ — '"+1." „^ ^"-"-i — '"."-1 



Now we have ?r " +1 '" = — — — and similarly 7r "'" 1 — , 

7Tn+l,n+l I n+1 „ +1 I n , n ln,n J 7T„, n 

Finally we use (|2~2]) . □ 
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Remark 5. Notice that if a is a symmetric measure, then its moments of odd order 
are equal to zero consequently following \2. 1)) and 12. 6\) we deduce that coefficients 
b n are all equal to zero. On the other hand coefficients a n can be expressed as 
functions of determinants of moment matrix as it follows from h2.6\) . On its hand 
coefficients a n determine completely orthogonal polynomials p%(x), i — 1, . . . ,n by 
il.9\) . i 1.1 0\) and {l.J$ and consequently matrices Lj(a), i = l,...,n which by 
its side determine moment matrix M„ . In other words we have an algorithm for 
regaining moments from the sequence of determinants of the major submatrices of 
the moment matrix. 

We have immediate observations concerning system of equations (|1.9[) - (|1.10[) 
which relate coefficients of 3-term recurrence to coefficients ir n i and consequently 
to the moments. 

Proposition 3. Let us denote rj n i — ix n i IlJ=i a j- i) ^ n > : a„ > 0. 

H) Vo,o = 1.1 n , n = 1. f° r n>0, 



Hi) Vn,n-1 = ~E"=0 h h 



2 



*») Vn,n-2 = Eo<,<j<n-l b i h j ~ Ej=l a y 

Assume that Vn > : bi = 0, then 

J if n=2k-l 

V) ^° ~ { UU 4-x if n = 2k ' k ~ X > 2 ' ' • ■ ■ 

vi) 

( if l = 2k + l 

(2J) Vn-l = j (-l) fc E X<n<...< jk <n-1 nLl«L if l = M ' 

^ jm+l-im>2,m=l,...,fc-l 

for k = 1,2,..., n>2k. 

Proof. Is shifted to Section [3] □ 

Remark 6. Note that coefficients j are the power coefficients of the so called 

"monic" orthogonal polynomials i.e. orthogonal polynomials with leading coefficient 
equal to 1. 

Remark 7. Notice that (— l) k ^2 i<ji<— i Ilm=i a j m can a ^ so ^ e u;r ^ _ 

+ 1 — Jttt >2,m — 1, . . .,fc— 1 

ten as 

n-2k+l n-2fc+3 n-1 

(-D fc E 4 E a %- E <■ 

j'i=i i2=ji+2 jfc=ifc-i+2 

As a corollary we get also the following recursive formula expressing moments in 
terms of 3-term recurrence coefficients a n and b n . 

Proposition 4. i) rrij = — Efc=i Vj-i k-i m k, 

If we assume that all coefficients b n — n > 0, then we have simplified version 
of the previous statement: 

ii) m 2 k-i = k = 1,2, m 4 = <z?(a? + af), m 2k = (Ej=i a|)™2fe-2 - 

Ej=2 7 72fe-l,2fc-l-2j/ i 2fc-2ji fc > 3. 

Proof, i) We use the (jl.6p and (|2.4|) which leads to the identity Vj > 1 

E^.fc mfc = °- 

fc=0 
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Consequently rrij = — J2k=o 1 lj,k m k- Now we utilize (|3.2[) and get: 

m j = - :// / i • - a2 j-iVj-2,k + Vj-i, k -i) m k 

k=0 

J-1 i-2 J-l 

J-l 
fc=l 

ii) By i) we have m 2 j = - EfcTi 1 V2j-i,k-i m k = ~ E£=i %j-i,2n-i m 2n since 
with odd j are equal to zero. Now we recall that f?2i-i 2j-3 = — Y^iLi 2 a i by 
Proposition [3l iv). □ 

2.1. Connection coefficients and Radon Nikodym derivatives. In this sub- 
section we will express the called connection coefficients between two sets of TV— orthogonal 
polynomials. So let us assume that we have two moment matrices Mjy (ex) and 
Mat (S) . Let Ljy (cx) and Ljv(<5) be their Cholesky decomposition matrices and 
{Pat (a;, a)} and {Pn(x,5)} respective sets of N— orthogonal polynomials. Then 
we have 

Lemma 1. We have 

P N (x,S) = L N 1 (S)L N (a)P N (x,a) . 
Proof. This formula follows simple observation that 

X„ = L n (a)P n (x,a). 
Then we apply Proposition [5] i). □ 

Following slight modification (ratio of densities is substituted by the Radon- 
Nikodym derivative of respective measures) of Proposition 1 iii) of [9] we deduce 
the following general statement concerning : 

Corollary 2. If ^(x) = l/Q r (x) where Q r is a polynomial of order r (positive on 
supp 5) then for N > r + 1 the symmetric matrix 

L N 1 (a)M N (S) (h-\a)) T 

is a V— ribbon' matrix i.e. its (i,j) entries such that \i — j\ > r are zeros. 

Proof. By the above mentioned Proposition we deduce that the lower triangular ma- 
trix L^ v 1 (a)L]v((5) is a V— ribbon' matrix. Then we have L A 7 1 (a)LAr(<5)(L^ 1 (cx)~Ln(o~)) T 

= L A 7 1 (a)MAr(<5) (L jv 1 (q!)) T . Then we use the fact that AA T is a V— ribbon' ma- 
trix iff A is a lower triangular V— ribbon' matrix. □ 

As a more interesting consequence of Lemma [T] we have important expansion of 
Radon-Nikodym derivative of two measures a « S. 

Theorem 1. Let the two measures a and S both having all moments be such that 
[^(x)) 2 dS(x) < co, where ( 



a << S and J (^(x)) 2 dS(x) < oo, where jjj(x) denotes their Radon-Nikodym 
derivative. Then 

dd 



, OO 

(2.8) -^(x)=J2E«P j (X,6)p j (x,5), 
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in L2(supp 6, J 7 , dd), where T denotes Borel sigma field o/supp<5. In particular 
(2.9) [Ax)) 2 d5(x)='£(E a p j (X,5)) 2 . 

Additionally when ^2j- >0 (E a pj(X, 6)) 2 ln(j + l) 2 < oo, we have S almost everywhere 
convergence. 

Proof. Let us notice that if we write p n (x, S) — Yh=o In ii^i a )Pi{ x i a ) then j n (8, a) 
= £j=o 4j ($) m j( a ) = I Pn{x,5)da{x) = E a p n {X,8) since we have (|2.4|l . Now 
assume that a « S and more over that Radon-Nikodym derivative ^g-(x) is square 
integrable with respect to measure S i.e. J (^(x)) 2 dS(x) < oo, then following [9] we 
see that ()2.8[) is true. Besides (|2.9[) follows Bessel equality of orthogonal series. If 
J2j>o(EaPj{X, S)) 2 ln(j + l) 2 < oo then we apply Rademacher-Menshov Theorem 
and get almost everywhere convergence. □ 



Example 1. As a corollary we will get famous Poisson-Mehler expansion formula 
( $2.1(J\) , below). In order not to repeat too many known details we refer the reader 
to [S] , |10] as far as the ideas and calculations are concerned and to [3] in order to 
get more properties of the mentioned below families of orthogonal polynomials. 

Namely we will consider the so called q—Hermite polynomials defined for \q\ < 1 
as H n (x\q) / \/[n] q \, where H n (x\q) are monic polynomials satisfying 3-term recur- 
rence given by (2.3) of [TO] . 

We used here traditional notation common in the so called q— series theory: [n] q 
= (1 - q n )/(l - q), for \q\ < 1 and [n]i - n, [n] q \ = \[ n j=1 \j] g , with [0],! = 1 (o)„ 

= Jl"=o 1 (^ ~ a 9 I ); ( so called q—Pochhammer symbol). 

One can consider also the case q — 1 obtaining simmilar results but for the sake 
of simplicity let us consider only the case \q\ < 1. 

It is known that q—Hermite polynomials are orthogonal for \q\ < 1, x G S (q) 
— {x € R : \x\ < 2/y/l — q} with respect to measure with the density fN(x\q) 
whose exact formula is not very important and which is given e.g. in [10) (formula 
(2.10)). Measure with the density fN{x\q) is our measure 5. It is also known (see 
same references) that the measure with the density : 

°° (1 - p 2 q k ) 
fcN (x\y,p,q) = f N (x\q) Y[ 



, _ Q w k (x,y\p,q)' 
where 

w k (x, y\p, q) = (l- p 2 q 2k ) 2 - (1 - q)pq k (l + p 2 q 2k )xy + (1 - q)p 2 (x 2 + y 2 )q 2k , 

for x,y € S(q), \p\ < 1 for \q\ < 1 has orthonormal polynomials equal to the so 
called Al-Salam-Chihara polynomials P n (x\y,p,q) satisfying the following 3-term 
recurrence given by formula (2.6) of [10] divided by {p 2 ) n [n] q \ as it follows from 
Proposition l,iii) of |10j (to get orthonormality) . 

Measure with density fcN it is our measure a. Following formula (4-7) in [?] we 
deduce that 

E a H n (X\q) = p n H n (y\q). 
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y G S(q) and \p\ < 1 are some parameters. Details are in [10] but can be traced to 
earlier works of Bryc, Matysiak, Szablowski. 

da (1 -p 2 q k ) 

Notice also that this function is bounded from above and as such square integrable 
with respect to any finite measure on S(q). Again details of the proof of this simple 
fact are in [10] . Now following K2. 8\) we get: 

00 / 1 2 k \ j 

,2 - io) E^^=z£< H ' mM - 

for every y £ S (q) and almost all x <E S (q) . Notice that for q = 1 12. 1 0\) is also 
true but it requires some more properties of Hermite polynomials. 

Remark 8. Situation described above is an illustration of the situation often met 
in the theory of Markov processes. Namely suppose that we have process X = 
{X t : t € T}, where T is some ordered set of infinite cardinality and Vi € T : 
X t is a random variable with support of infinite cardinality. Suppose dP t is the 
distribution of X t and that E t \X t \ n is finite for all t and n. Suppose also that |pn "* | 
are polynomials orthogonal with respect to dPt . Further suppose that the conditional 
distribution of X s given Xt — y for s > t i.e. dC s ,t is absolutely continuous with 
respect to dP s and that (x) is square integrable with respect to dP 8 for every 
s > t and y £ suppXt- Then as it follows from Theorem^ in L2 (supp X s , J 7 , dP s ) 
we have: 

^{x)dP s = (Y,E s Jf ) (X s )p^\x))dP s . 

That is we get expansion of the transfer function of our process. 

2.2. Linearization coefficients. Notice that Propositions Q] and [5] allow us to 
formulate an algorithm to get so called 'linearization coefficients'. Let us recall 
that linearization formula is popular name for the expansions of the form 

Pn (^) Pm (^) — ^ ^ Cn,m, (^) • 

i=o 

The problem is to find coefficients c n _ m j for all n, m > 1. One can easily deduce 
(following general properties of orthogonal polynomials that c njTnj = for j = 
0,...,|m — n\ — 1. Let us fix N > and let us consider vectors and matrices as 
introduced in the Introduction. We have 

XatX^ = AArPArP^-A^-. 

Now notice that PjvP^ is the (N + 1) x (N + 1) symmetric matrix with (i,j) 
-entries equal to pi (x)pj (x) , while XjyX^ is the (N + 1) x (N + 1) symmetric 
matrix with -entries equal to Xi+j. Following (ll.5[) we know that: 

i+j 
k=0 
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Hence introducing symmetric tensor T2N with entries A»+j,fc, i,j — 0,...,N, k 

SV 



0, . . . , N we have XjvX^ = T2jvP2A'- So we deduce that 



P^Pw — njvT 2 Arn^P 2 Ar. 



Another words 



(x)p rn (x) = y y Kn^TTmjX. 

k 0<?<n, 



/ 



Following general properties of orthogonal polynomials we deduce that Vfc < |n • 



'" : I S 0<i<n, TTn^TTmjAi+j^ ) = 0. 
0<j<m,j+i>k 



3. Proofs 

Proof of Proposition [H i) Follows uniqueness of both Cholesky decomposition and 
orthonormal polynomials provided sign of the leading coefficient is selected, ii) We 

J P n {x,a)Pn{x,a)da(x) = L" 1 J X n X^da(x) (L^ 1 ) 7 " 

= L-^n {L- 1 ) T = In- 

ii) We have 

P^s, a)P„(y, a) = (L- 1 ) T L" 1 Y„ = (X„) T (L Tl I^)- 1 Y B . 
We obviously have 



X„r/A„ < X^M-^n < |X n |VAi, and |X„| 



i=0 



in) 



Pl(x)M n P n {x)da 



tr{M n P n {x)Pl{x))da 



= trM„L" 1 M„ (L,; 1 ) 7 = trM n . 



iv) We have by Proposition^ i) £ P n ( e *)pr( e -«)dt = JL / o 27r e^M^e^- 
where we denoted e^(i) = (1, e lt , . . . , e mt ). Secondly notice that — th entry of 



the matrix P„(e rt )P„ (e lt ) is equal XX=o 7r »,fc 7r i,fc e e which , by the same 
assertion is equal to (i, j)— th entry of the matrix M^ 1 . Second statement follows 

the fact that £ £™ =0 JcT |Pj( eit )| 2 d * is the trace of h lT Pn(e u )P^(e- lt )dt. 

v) By Proposition [21 ii) considered for x — y = O.We get X^ob'WI 2 = 
O^M-^n, where 0^ = (1, 0, ... , 0), which means that £" =0 b*(°)| 2 is (°> °) entr y 
ofM" 1 . 

vi) Following (jl.lll) we see that q n (0) = Sj=i ^n,jfnj-i- Another words Q„(0) 

= Tl n fi n . Now 5Z™ =1 ?j (0)) 2 = Q^(0)Q n (0). On the way we utilize the fact 

that n^n„ = M" 1 . Following similar arguments we have j=i Qj (0) = 
(l,0,...,0)M„/x„. 
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vii) Let us denote p n (x, a) = ^ n+1 i 1 g ^(„ +2 ) Y^i=oPi-i x i a )- ^ satisfies recursion 



/ n + 1 log (n + 2) , — — 2 

p„ + i(x,a) = W — — - — J- — pn(x) + p n +i [x) Vn + 2 log (ra + 3). 

V n + 1 log (n + 3) 

Since we have 5~L-^ n ? — rrf /"^ 2 \ n < oo we deduce by Rademacher-Menshov the- 

^^n>U (n+1) log 4 (n+2) J 

orem that series J2n>o ^ n +i\o^(n+2) conver g es a— a.s. Further we apply [5] (Thm. 
5). X " " □ 

Proof of Proposition \1.7\ Multiplying both sides of (11.91) and (jl.101) by 11™= l a * we 
see that quantities rj satisfy the following system of equations: 

(3-1) Vn+1,0 + b nVn,0 + a l 7 ln-l,0 = °> 

(3-2) ^n+lj + bnVnJ + a n r ln-l,3 = Vn,j-li 

n > 0, j < n + 1. i) Follows immediately formulae (|2.6[) . 

ii) Take j = j + 1 in^£2§- Since 7? fc >n+1 = for fc < n + 1 we get V n +i, n +i = 1n,n 
with n 00 = 7r ,o = 1. 

iii) Again we consider (|3.2p with j = n which takes the following form: 

and notice that for n = we have 7Ti,o + &o = 0. Further we deduce by induction. 
Suppose iii) is true for j < k. Let us consider n = k + 1 we have »/ fc +i ^ + b k = 
— ^2j = q bj from which it immediately follows our assertion. 

iv) Notice that from p. II) considered for n = 1 wc get: rj 2 + 61?^ + af = so 
r/ 2 = b()bi — af so our formula is true for n = 2. Further we deduce by indiction. 
Assume that it is true for j < n—1. Considering (|3.2|) for j = n— 1 we get rj n+1 n _ 1 
+ b n rj n „_j +ct„^„_i „_i = Vnn-2- Now we put formulae from assertions ii) and iii) 
and induction assertion and get: r /n+l n _ 1 - b n Y^Zq bi + a 2 n = J2o<i<j< n -i b ibj ~ 
Ej=i a )- Now since Eo<«<j<n b i b i = b n E"=o b i + Eo<i<j<»-i &i&j,and X)"=i a| 
= Y^j=i a j + a n we S et our assertion. 

v) Let us consider (|3.1[) with 6„ = 0. We get then r] n+10 = ~af l r] n _ 10 . Recall 
that then rj = 1 and 77 10 — 0. So we see that rj n with odd n must be equal to 
zero. 

vi) To see that under our additional assumption rj n n -2k+i = 0' ^ = 1,2,..., 
n > 2/c — 1 is easy since then our formula (|3.2[) becomes now: 

(3-3) 7 ln+l,n+l-2k = — a n 7 7n-l,n-2fc+l + 7 ?n,n-2fc- 

Besides we have r] 2k+1 = by v) and t?„_j „_2/c — by induction assumption. 
Hence let us consider the case of even differences in indices i and j in 77^ ^ . 

The proof will be by induction. First notice that since sign of r) nn _ 2k is (— 
and of «„_i.„_2fc+i is (~ l) fe_1 by induction assumption we deduce that the sign 
°f Vn+i n+\-2k ^ s (~ -0 as claimed. Secondly notice that 7y n „„ 2 fe being a sum of 
certain products of k different a? does contain only products with j < n, hence 
the two sums i.e. — dnVn—i n-2k+i and r\ n n _ 2k do not contain similar products 
and moreover since i] n _ 1 „_ 2 fc+i does contain product of k — 1 different cn with 
j < n — 1. Consequently 77„+i i , i +i_2fe i s the sum of products of /c different en with 
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j < n+ 1. Secondly notice that the sum in (|2.7j) contains ( n k k ) summands. Hence 
number of summands in (|3.3|) is equal to 

f n - 1 - (k - 1)\ (n - k\ _ fn + 1 - fc N 
k- 1 Jv fc /V & 
by the well know property of the Pascal triangle. □ 
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